function [ s ] = evaluateQuadraturePoint( p, q, qn )
%EVALUATEQUADRATUREPOINT Summary of this function goes here
%   Detailed explanation goes here

r = p - q;

rn = norm(r);

rg = r / rn;

qnc = [qn(2); -qn(1)];

M = [rg(1)*rg(1), rg(1)*rg(2); rg(2)*rg(1), rg(2)*rg(2)];
N = [0, dot(rg, qnc); -dot(rg, qnc), 0];
I = diag(ones(2, 1));

prpn = dot(r, qn) / rn;

s = -1 / (4 * pi * rn) * (prpn * (I + 2 * M) - N);

end

